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ABSTRACT 

We show that the star-forming regions in high-redshift luminous and ultraluminous infrared galaxies 
(LIRGs and ULIRGs) and submillimeter galaxies (SMGs) have similar physical scales to those in 
local normal star- forming galaxies. To first order, their higher infrared (IR) luminosities result from 
higher luminosity surface density. We also find a good correlation between the IR luminosity and IR 
luminosity surface density in starburst galaxies across over five orders of magnitude of IR luminosity 
from local normal galaxies to z ~ 2 SMGs. The intensely star- forming regions of local ULIRGs are 
significantly smaller than those in their high-redshift counterparts and hence diverge significantly 
from this correlation, indicating that the ULIRGs found locally are a different population from the 
high-redshift ULIRGs and SMGs. Based on this relationship, we suggest that luminosity surface 
density should serve as a more accurate indicator for the IR emitting environment, and hence the 
observable properties, of star- forming galaxies than their IR luminosity. We demonstrate this approach 
by showing that ULIRGs at z ~ 1 and a lensed galaxy at z ~ 2.5 exhibit aromatic features agreeing 
with local LIRGs that are an order of magnitude less luminous, but have similar IR luminosity 
surface density. A consequence of this relationship is that the aromatic emission strength in star- 
forming galaxies will appear to increase at z > 1 for a given IR luminosity compared to their local 
counterparts. 

Subject headings: galaxies: evolution — galaxies: high-redshift — infrared: galaxies 



1. INTRODUCTION 

Luminous and Ultraluminous Infrared Galaxies 
(LIRGs and ULIRGs) whose total infrared luminosities, 
L(TIR), are in the range of 10 11 - 10 12 and > 10 12 L , 
respectively, are among the most important populations 
for the study of galaxy evolution. Their extreme L(TIR) 
signifies intense star formation hidden by dust and only 
visible through the infrared (IR) emission resulting 
from the reprocessing of UV photons that originate 
from populations of rapidly forming hot young stars. 
They are among the most extreme star-forming regions 
known. 

A majority of the local ULIRGs are disturbed systems 
of galaxies undergoing mergers. Some ULIRGs harbor 
active galactic nuclei (AGN) often thought to be asso- 
ciated with merger activity and the resulting enhanced 
dense molecular ga s fractions observed in the nuclei of 
these systems (e.g. iGao fc Solomon! l2QQ4t Uuneau et al.l 
2009) . Locally, lum inous IR galaxies are very rare 
(jLagache et al.l l2QQ5h but they contribute significantly 
to th e total IR energy density beyond redshift z ~ 1 
(e.g., iLeFjoc'h e t al. 2003 iPerez-Gonzalez et al.l 120051 : 
lR,odighiero et al.l l201Q[ ). At high redshift, z > 2, lu- 
minous IR galaxi es are typified b y submillimeter galax- 
ies (SMGs; e.g. iBlain et al.ll2002f ). So far it is common 
to view SMGs as a more luminous counterpart of local 
ULIRGs. This naturally leads to th e view of SMGs bein g 
disturbed, interacting systems (e.g. iTacconi et al. 2008). 

However, there are many indications that the SEDs 
of the high redshift infrared galaxies differ systemat- 
ically from those of local galaxies of similar luminos- 
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ity. iPapovich et~a71 (|2007D found that the 70 /im (ob- 
served) outputs tended to be weaker relative to those 
at 24 /im (observed) than expected from local tem- 
plates. Many authors have found that the aromatic 
bands in these galaxies appear to be similar in structure 
to th ose in significantly lower- luminosity l ocal galaxies 
fe^ lRigbv et al.ll2008t iFarrah et al]l2QQ8t iTakagi et all 
2010). The far infrared SEDs appear to be cold, again 
similar to those o f lower luminosity local galaxies (e.g.. 
Pope et all 120061 : iSvmeonidis et al.l 120091 iMuzzin et all 



2010f ). These findings suggest an underlying physical 
difference between local l uminous infrared g alaxies and 
those at high redshift. Rig bv et"aTI (j2008h suggested 
that such a difference might arise either through reduced 
metallicity or lower optical depth due to a g reater ex- 
tent of the emitting regions. lErb et ahl ([2006) find that 
the metallicities are of o rder three times lower at z ~ 2 
for a given galaxy mass. lEngelbracht et all (|2008f ) show 
that, for local galaxies, there is only a weak trend in 8 
/im vs. L(TIR) down to 1/3 solar metallicity, and that 
at lower metallicity the 8 jam luminosity is suppressed. 
This correlation is also reported for z ~ 2 galaxies by 
iReddv et al.l (j2010f ). We conclude that reduced metallic- 
ity is unlikely to be the primary cause of the changes in 
SED with redshift. We therefore turn our attention to 
the second possibility, that the galaxies have structures 
different from local ones of similar luminosity. 

Recent high-resolution studies of SMGs in the sub- 
millimeter, radio, and near-IR have shown that their 
star-forming regions are generall y relatively extended , 
with diameters of or der 1 — 10 kpc ( Chapman et al. 2004 
Mux low et aTll2005l: iBiggs fc Ivison 2008; Both well et"at | 
20101 : ICasev et al . 2009; Iono et al. 2009; Leh nert et all 
2009; CarillietaD 120101 : Swinbank et al. I2010at 
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iTacconi et all 120101 : lYounger et all 120101 ) . Thus, 
the surface densities of the star formation rate (SFR) 
must be substantially lower in the high-z galaxies than 
for the local ones of similar L(TIR). This paper explores 
the importance of this difference in explaining the 
different SED behavior. We make the simplest possible 
assumption, that the SED is a function of the SFR 
surface density, Ssfr, and that galaxies with similar 
Ssfr, and consequently similar L(TIR) surface density, 
Sl(tir)? will nave similar optical depths in their star 
forming regions and similar SEDs. In section [2j we 
describe the compilation of data for the study as well as 
discuss our derivation of a consistent set of sizes for the 
star forming regions for < z < 2.5. Section [3] shows the 
results and discusses the implications for both galaxy 
evolution and for improved estimates of the SFRs of high 
redshift IR galaxies. We assume a ACDM cosmology 
with n m 0.3, Q A 0.7, and H = 70 km s^Mpc" 1 
throughout this paper. 

2. THE MEASUREMENTS 

A major challenge in studying the Xl(tir) of galaxies 
across a broad redshift range is to obtain physical size 
information for the IR-emitting regions in a consistent 
way for both the local and high-redshift samples. We 
discuss in this section the compilation of physical sizes 
and IR luminosity data from the literature and our steps 
to compare them on the same metric. Despite the abun- 
dantly available high-quality data for local galaxies, it 
is necessary to select a subset that can best match the 
high-redshift sample. We therefore will begin our discus- 
sion with the high-redshift compilation, then follow with 
the local compilation. 

2.1. Intermediate and High- Redshift Compilation 

The past five years have seen an unprecedented num- 
ber of high-resolution, multiwavelengths observations of 
high-redshift galaxies. Our compilation of the physical 
sizes of intermediate and high-redshift starburst galax- 
ies is based on radio interferometric observations at 
the Multi-Element Radio Linked Interferometer Network 
(MERLIN) and the Very Large Array (VLA), as well as 
submillimeter interferometric observations at the IRAM 
Plateau de Bure Interferometer (PdBI). Submillimeter 
and 1.4 GHz observations, rather than rest-frame op- 
tical or near-IR, are chosen because they closely trace 
the star-forming regions of galaxies while being little af- 
fected by the old stellar light. The submillimeter is a 
more direct tracer in this regard because it observes ther- 
mal emission from dust while the 1.4 GHz data probe 
synchrotron radiation from super novae remnants, whose 
progenitors are short-lived massive stars and hence indi- 
rectly trace star-formation. B ecause of the radio- i nfrared 
relation within galaxies (e.g.. iMurphy et a"T1l2006l ). radio 
observations should give a valid measure o f the size of 
the act ive region in a star forming galaxy. Ilvison et al.l 
(|2010al ) report that the IR-radio relation could evolve 
at high redshift s. However, their reported evolution law 
given by (1 + ^)-° 15±0 - 3 i s small and unlikely to affect 
significantly our estima tion of IR-emitting re gion size us- 
ing radio observations. iSargent et al.l (J2010f). al o ng wit h 
a recent result using Herschel bv Ilvison et al.l (|2010bl ). 
also found little or no evolution in the IR-radio relation 
at z < 2. 



Although submillimeter transitions (e.g. CO) are more 
direct tracers of molecular clouds and star-forming re- 
gions than is radio emission, they tend to trace colder 
gas. In comparison, the star-forming regions in most 
local LIRGs and ULIRGs are sub-kpc in size (e.g., 
iCondon et al.lll99lh . Thus the physical size given by sub- 
millimeter observation tends to be systematically larger 
than that of radio observations. Moreover, selection of 
high-redshift submillimeter galaxies will tend to result in 
a sample with colder SEDs, which are brighter in submil- 
limeter wavelengths, and hence an inherently physically 
extended sample. This bias towards larger physical size 
should be more pronounced for the low— J CO transi- 
tions such as CO (1 — 0) and CO (2 — 1) and thus we 
will adopt sizes from higher— J CO transitions such as 
CO (3 — 2) and CO (6 — 5) if available. To be conserva- 
tive, we will take the radio sizes to be the physical size of 
the star-forming region in galaxies and treat submillime- 
ter sizes as upper limits to isolate this possible selection 
effect. 

2.1.1. High-z Submillimeter Compilation 

The subm illi meter data were c o mpiled from 
Tacconi et all (12006ft . ITacconi et all ([2010ft . and 
Dad di et al.Pj2010ft . These authors used the IRAM 
PdBI to study submillimeter-selected sam ples using CO 
transi tions with 0'/6 to ~ l'/0 resolution. ITacconi et al.l 
(200 61) observed s ix SMG s, including fou r sources from 
IGreve et al.l (|2003ft and iNeri etHI (|2003f ). and two new 
sources at redshifts 2.2 < z < 3.4. Their observations 
yield four resolved SMGs with a median FWHM diame- 
ter of ~ 4 kpc. One of the ITacconi et al.l ([20061 ) sources 
needs to be excluded because its radio emission is signif- 
icantly stronger than predicted by the IR-radio relation , 
indicating the presence of an AGN. ITacconi et all (|2010ft 
reported the physical sizes (effective CO diameters) of 
two sources in the Extended Groth Strip at redshifts 
z ~ 1.1 to be 13 and 16 kpc in diameter, which they 
have found to be smaller but consistent with the size 
based on an I-band observation. It should be noted 
that although ITacconi et al.l (|2006f ) find the SMGs to be 
"compact" , they still generally are a factor of 100 larger 
in area than local ULIRGs of s imilar IR luminosities. 
Additionally, iDaddi et al.l (|2010ft reported sizes of four 
BzK-selected z ~ 1.5 star- forming galaxies (one member 
of this samp l e also has 1.4 GHz size measurement from 
iCasev et al.l ([2009ft , which we adopt in preference to the 
CO data). In total we have eight data p oints from the 
submi llimeter ob servations; th r ee fro m ITacconi et al.l 
(2006 |), two f rom ITacconi et al.l (|2010), and three from 
IDaddi et al.l (|2010ft . IR luminosi ties for thes e eight 
galaxies were estimated via the Rieke et al. (2009) 
forma lis m based on 24 /im fluxes from lHainline et al.l 
([2009ft . iMenendez-Delmestre eFall ([2009ft . and Spitzer 
deep imaging in Great Observatories Origins Deep 
Survey (GOODS, Dickinson et al., in prep.). 

2.1.2. High-z Radio Compilation 

The 1.4 GHz ra di o ob se rvations were draw n 
from IChapman et al.l (|2004f ): iMuxlow et "all p005); 
Biggs fc Ivisonl (|2008ft and ICasev et al.l (|2009). 
ChapmaneFaD (|2004f ) used MERLIN+VLA with 
0'/3 resolution to study 12 SMGs selected from the 
HDF at a median redshift of z = 2.2 ± 0.2. 8/12 of 
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Fig. 1. — To determine the physical size of star-forming regions 
consistently for both local and high-redshift galaxies, we need to 
degrade the image of local galaxies to resolution similar to that ob- 
served at high-redshift. We convolve local images with a Gaussian 
beam of 4 kpc at the galaxy's distance and fit the convolved images 
to estimate the deconvolved diameters. A typical case, NGC 3627, 
is shown here with the original Spitzer 24 /im image at the left panel 
of the inset. Subsequently we show the convolved images with the 
corresponding sizes of the convolved Gaussian beam and the 2D 
Gaussian fits to the convolved images overplotted as ellipses. This 
example demonstrates that the estimated deconvolved size is rel- 
atively insensitive to the choice of Gaussian used to convolve the 
original image. 

their targets were resolved with an effective median 
diameter measured ab ove the 3a contour of 7.0 ±1.0 kpc. 
IBiggs fc Ivisonl J2008) used MERLIN+VLA to observe 
12 SMGs in the Lockman Hole atl.2<z<2.7 with 
0'/2 — 0'/5 resolution (all sources were resolved). The 
sample was selected from SMGs foun d by surveys car- 
ried out with SCUBA an d MAMBO (Scott e tld][200l 
iGreve et al.l 120041 : iCoppin et al.l 120061 ). IBiggs fc Ivisonl 
(2008) reported physical sizes of these SMGs ranging 
from 1 — 8 kpc with a median of 5 kpc. Lastly, 
ICasev et al.l (|2009h used MERLIN (0'/3 - 0'/5 resolution) 
to observe seven ULIRGs and a HyLIRG selected by 
optical and radio color criteria from GOODS and the 
Lockman Hole at redshifts 0.9 < z < 2.4. They resolved 
all eight and report an average physical diameter of 
5.0 ± 1.2 kpc, corresponding to the surface area within 
t he 3a detection level. 

iMuxlow et all ([20051 ) used MERLIN+VLA to conduct 
a deep 1.4 GHz survey of the HDF and Hubble Flank- 
ing Field (HFF) and studied 92 radio sources at 0.18 < 
z < 4 with 0'/2 — 0'/5 resolution. All but one source 
was resolved. For t his study, we exclude sources that 
IMuxlow et all (|2005l ) classified as AGN or AGN candi- 
dates. We also exclude the Muxlow sources that are 
observed only by the VLA and hence have lower reso- 
lution that may affect the accuracy of size estimates, as 
well as those with complex morphology s uch that their 
size e stimation req uires visual inspection (jMuxlow et al.l 
20051: iThrall et al.l l2QQ7h . Out of 92 objects from the 



Muxlow e t al. (|2005f ) sample, 72 have spectroscopic red- 



shifts and 27 pass the aforementio ned criterion and have 
a 24 /im counterpart in GOODS. IMuxlow et "all (j2005D 
reported sizes in terms of t he largest angular s ize deter- 
mined by Gaussian fitting (jThrall et al.l 12007). Taking 



the largest angular size, which is the size of the major 
axis, as the diameter directly would overestimate the sur- 
face area; we assume that the sizes of the minor axes of 
these galaxies are 0.6 of the major-axis (the largest an- 
gular FWHM) and then calc ulate the surface area s and 
circularized diameters for the IMuxlow et al.l (j2QQ5h sam- 
ple as reported in Table El This value is an average ratio 
of min or-to-mai or axes of SMGs from the Bigg s fe Ivisonl 
(200 81) sample . An average physical diameter for the 
IMuxlow et all (120051) sample is 4.5 kpc. L(T I R) fo r 
the [Chapma net al.l (12004J) . the IMuxlow et all $2005), 
and the ICasev eHa] (l2009|) samples in the G OODS field 
were estimated via the [Rieke et all (|2009[ ) formalism 
and the Spitzer 24 a m im aging in GOODS. L(T I R) for 
the IBiggs fc Ivisonl ([2008) and the ICasev et al.l (|2009D 
samples in other fields are also estimated using the 
iRieke et"aTI (120091) formalism b ased o n the 24 am fluxe s 
reported bv llvison et al.l (j2007l ) and ICasev et al.l (j2009f ). 

Apparently there are two approaches in reporting 
galaxy sizes from the radio measurements: the de- 
convo lved FWHM (jMuxlow et al.l 120051 : IBiggs fc Ivisonl 
2008) and the circularized diameter en closing the sur- 
face area of the > 3a r adio detections ([Chapman et al.l 
120041 : ICasev et al.ll2009D . Despite the methodological dif- 
ferences between these two approaches, they agree very 
well for the objects with overlapping observations, espe- 
cially in the HDF where there are six such galaxies. This 
agreement is expected, since the area detected above 3a 
is likely to be consistent with that of half-max power 
given the typical levels of signal to noise ratio in the ob- 
servations. We also note a very good agreement of sizes 
from 1.4 GHz and CO (3 - 2) observations of a HDF 
object, SMMJ 1 23707 +62 14SW, where o b serva tions from 
IMuxlow et all (|2005l ) and lTacconi et al.l (|2006f ) overlap. 

Although galaxies drawn from the aforementioned 
samples are likely dominated by star-forming activity, 
we independently confirm this by inspecting whether 
their ratio of 850 jam and 1.4 GHz fluxes follows 
the IR-radi o rela tion for starburst galaxies given by 
Riek e et"a"T1 (j2009[ ). Galaxies with radio flux signifi- 
cantly stronger than predicted by the IR-radio rela- 
tion likely harbor radio-loud AGN. Indeed, we have 
found that the only two objects with F850/F1.4 ^> 
2 are the most compact objects in our compila- 
tion (consistent with AGN-domination) and thus we 
exclude these objects (SMMJ163650.43+405734.5 and 
SMMJ105207.49+571904. q). It is worth n oting that this 
test is unnecessary for the IMuxlow et al.l (j2005f ) sample, 
where both radio spectral slope and radio morphology 
are already employed to identify AGN and AGN candi- 
dates. 

In total our compilation has 48 intermediate and high- 
redshift starburst galaxies with a median redshift of 1.0 
and a median size of 5.1 kp c. Our primary samp le com- 
prises the 27 galaxies from IMuxlow et al.l (j2QQ5h . which 
should provide an unbiased and complete sample, se- 
lected by radio flux alone. The results from the ad- 
ditional 21 ga l axies are consistent with those from the 
IMuxlow et al.l (|2005f ) observations. 

2.2. Local Compilation 

For local galaxies, there is no homogeneous set of radio 
or submillimeter images suitable for our needs. There- 
fore, we use two additional measures of the star forma- 
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tion rate, Paschen-a images (|Alonso-Herrero et al.l l2006) 
and images at 24 /am from a number of sources. As with 
the radio, both measures are not strongly affected by 
extinctio n and are good tracers of star forming activ- 
ity (e.g., iMuri^ilaD HQS!; [CaliettreF^ (S^Oj) . Al- 
though radio images of local galaxies have somewhat 
larger extent than these other indicators, the high surface 
brightness areas that would dominate the high-redshift 
measurement have v ery similar morphologies and sizes 
(jMurphy et al-ll2QQ6h , 

The high redshift galaxies have limited structural infor- 
mation, typically equivalent to a FWHM from Gaussian 
fitting or deconvolution of the observed image, in which 
the galaxy is only modestly well resolved. The FWHM 
of the beam is typically 0'/3 to l'/O. A l'/O beam cor- 
responds to a diameter of 5 kpc at z = 0.4, 7.5 kpc at 
z = 0.8, and 8.5 kpc at z = 2. However, the images of lo- 
cal galaxies often provide many more resolution elements 
across the target. To put them on the same scale as the 
radio images at high redshift, we convolved them with a 
Gaussian beam, and then determined sizes by Gaussian 
fitting and deconvolution. We demonstrate with a typi- 
cal case of local star-forming galaxy in Fig. [T] that this 
method yields a robust measurement of physical sizes for 
local galaxies. 

Despite our efforts to put all the images on the same 
basis, a possible source of systematic error remains. An 
image with low signal to noise will tend not to go as 
far into the wings of the star-forming activity and hence 
there will be a tendency for a systematic reduction in the 
estimated FWHM with decreasing signal to noise. Thus, 
the sizes of the high-redshift galaxies may be underesti- 
mated compared with the local ones. Since the major 
result of this paper is that the high-z galaxies have sub- 
stantially larger sizes than local ones of similar luminos- 
ity, the result is that we may understate this conclusion, 
not that it would be undermined. 

We describe our compilation of local IR galaxies in 
detail in the Appendix. It comprises 19 normal star- 
forming galaxies (L(TIR) < 10 11 L©), 21 luminous IR 
galaxies (LIRG, 10 11 < L(TIR) < 10 12 L ), and 4 ultra- 
luminous IR galaxies (ULIRG, L(TIR) > 10 12 L ). 

2.3. HST Imaging of the High- Redshift Sample in the 
GOODS-N Field 

Since 36 of 48 high-redshift galaxies in our sample are 
in the Hubble Deep Field imaged by the Great Observa- 
tories Origins Deep Survey (GOODS; Dickinson et al., in 
prep.) using the Advanced Camera for Survey (ACS) on 
the Hubble Space Telescope, we investigate the GOODS 
imaging (Version 2.0) of these galaxies in the B, V, z, and 
z bands. Details about the ACS observa tions of GOODS 
as we ll as the data can be found in Giaval isco et al.l 
(2004) and at the survey's websit^E 

GOODS optical images are used to confirm indepen- 
dently that the sizes and positions of the star-forming 
regions determined by radio or submillimeter observa- 
tions are consistent with the high-resolution optical ob- 
servations. The synthesized beam sizes of the radio and 
submillimeter observations are Of! 2 — Of! 5 while the optical 
imaging in GOODS/ ACS has resolution at ~ O'.'05/pixel. 

1 http:/ /www. stsci.edu/science/goods 



Although the radio and submillimeter observations can 
constrain the size of the star formation regions without 
being affected by stellar emission, their interpretation 
benefits from high-resolution optical imaging, especially 
for morphological classifications. We illustrate the circu- 
larized size for the star-forming regions compared to the 
optical extents of their host galaxies as seen in various 
optical filters in Fig. O The star-forming regions seen 
at radio and submillimeter wavelengths typically coin- 
cide with the central part of the optical structure, but 
there are a few cases where the star-forming region is 
at the collision interface of an interacting system (e.g. 
J123714+621558, J123716+621643). We will discuss our 
qualitative assessment of optical morphologies to address 
the structure of high-redshift star-forming galaxies in 

3. RESULTS 

Our determinations of the physical sizes and luminosi- 
ties for star-forming galaxies from the local sample out 
to the high-redshift SMGs are summarized in Fig. 03 
We find the physical sizes of star-forming galaxies to be 
comparable within an order of magnitude across the en- 
tire IR luminosity (10 9 — 10 14 L©) and redshift range 
(0 < z < 2.5). Local LIRGs and ULIRGs, however, are 
significantly smaller in size. It follows from Fig. [3] and 
is apparent from the diagram of Sl(tir) as a function of 
L(TIR) in Fig. |4] that there is a correlation extending 
more than five orders of magnitude between Sl(tir) and 
L(TIR). A linear fit to this correlation yields a formal 
slope of 0.9. However, given the inhomogeneity of the 
methods employed to estimate uncertainties in each of 
the subsamples comprising our compilation, we expect 
the uncertainty for this formal slope to encompass the 
slope of unity and hence suggest that the two parame- 
ters are proportional to each other. 

Local LIRGs and ULIRGs have 1 — 4 orders of mag- 
nitude higher Ul(tir) than indicated by this correlation 
suggesting that the starburst regions in these galaxies 
are not representative of their high-redshift counterparts 
of similar L(TIR). In fact, local ULIRGs with L(TIR) 
~ 10 12 have a comparable Sl(tir) to high-redshift SMGs 
with L(TIR) - 10 14 L . 

3.1. Ssfr as a Tracer of Star- Forming Environment 

The aforementioned correlation suggests an explana- 
tion for the strong aromatic emissions as well as the 
SEDs of star- forming galaxies at high redshift. If we 
assume that Ssfr, and hence El(tir)? * s an accurate 
indicator of the star-forming environment in starburst 
galaxies, Fig. |4] would suggest that z ~ 1 ULIRGs have 
similar environments to local LIRGs with an order of 
magnitude lower L(TIR), and likewise z ~ 2 SMGs have 
similar environments to local ULIRGs with two orders of 
magnitude lower L(TIR). We test this possibility with 
three independent approaches. 

First, we compare the logarithmic ratio of 24 /xm flux 
and 1.4 GHz flux, log(S24/Si.4) or #24, predicted by the 
SEDs of local LIRGs and ULIRGs and the observed flux 
ratio for high-redshift galaxies. At redshift 1 < z < 2.5 
the Spitzer 24 /im observed bandpass probes rest-frame 
wavelengths of 7 — 12 /im and thus <?24 at these redshifts 
is highly sensitive to emission from aromatic complexes, 
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Fig. 2. — HST ACS imaging of the subsample of our compilation in the GOODS field in B, V, i, and z. The circles indicate the size of 
the circularized diameter for the star-forming region in the corresponding galaxy. The bar at the lower left corner of the images in the right 
column represents a physical scale of 5 kpc. We found that ~ 2/3 of the subsample in GOODS display isolated, quiescent star-forming 
galaxies while a few cases show obvious signs of galaxy interaction as well as disturbed morphologies that can be associated with asymmetric 
clumps of star formation, especially at high-redshift. [See the electronic edition of the Journal for Fig. 02-0^.] 



which are sensitive to the star formation environment. 
We first investigate the intrinsic scatter of #24 by consid- 
ering the scatter observed locally derived from the IRAS 
and VLA observations, shown in the top panel of Fig. 
[U If we assume a similar scatter around the mean value 
of #24 and project this scatter to higher redshifts, the 
bottom panel of Fig. [5] demonstrates that the observed 
scatter at high redshift is within the expected scatter 
based on local observations. 

We then show in the right panel of Fig. [5] that #24 
for high-redshift galaxies is consistent with the ratios 



from local LIRGs and ULIRGs with significantly lower 
L(TIR). Consider the SMGs at z - 2.4 with L(TIR) 
~ 10 13,5 L ; the correlation in Fig. |4] suggests their 
starburst environment to be similar to local ULIRGs with 
L(TIR) ~ 10 12 Lq and #24 is clearly inconsistent with the 
extrapolation for local ULIRGs at L(TIR) - 10 13 5 L© 
while agreeing with the SED for L(TIR) - 10 12 L lo- 
cal ULIRGs. Another way to interpret this is illustrated 
by Fig. [6] that the aromatic features are significantly 
stronger compared to local SED templates at redshifts 
greater than z ~ 1. 
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Fig. 5. — (Top) The logarithmic ratio of 24 /im and 1.4 GHz 
flux densities, ^24, for high-redshift star-forming galaxies {right 
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The local distributio n is from the IRAS Bright Gala xy Catalog 
(SanderseFaD[2003) and the NRAO VLA survey (| Condon et al.l 
1998). The thick centr al line in the righ t panel shows the track 
of (?24 predicted by the Rieke et al. (2009) SED for a galaxy with 
L(TIR) of 10 12 Lq, with dark and light grey shades showing the 
la and 2a extents assuming the local distribution of ^24- Despite 
a relatively large scatter, the distribution of q24 at intermediate 
and high redshifts broadly agrees with the scatter observed lo- 
cally. (Bottom) The l ogarithmic ratio of 24 fim flux and 1.4 GHz 
flux predicted by the IRieke et all {2009) local LIRG and ULIRG 
SED templates (solid lines) compared to the observed ratios at 
0-5 < z < 2.5. The observed ratios are consistent with the predic- 
tions from local templates for galaxies with significantly lower IR 
luminosities. For instance, galaxies with L(TIR) of 10 13 — 10 14 Lq 
(blue dots) at z ~ 2.5 have a ratio consistent with the template for 
L(TIR) ~ 10 12 L Q (green line). 



Second, we consider a result from lRigbv et al.l ((2008), 
specifically their Fig. 3, that shows aromatic emissions 
for SMM J1635554.2+661225, a lensed SMG at z = 
2.516. This SMG has L(TIR) of 10 119 L , which would 
be hardly detectable if not for the 22 x lensing magnifica- 
tion. The magnification also allows high signa l-to-noise 
mid-IR spectroscopy with Spitzer, in which lRigbv et ahl 
(j2008[ ) found aromatic emission band shapes similar to 
those of NGC 2798 and to the average mid-IR spectra of 
13 local nor mal starburst galaxies with average L (TIR) 
of 10 10 - 7 L Q dBrandl et all 12003 iDale et aLll2007h . This 
result is consistent with the relationship given in Fig. |4j 



high-redshift galaxies with L(TIR) of ~ 10 119 L would 
have Sl(tir) similar to local galaxies with an order of 
magnitude lower L(TIR). 

Third, we compare a stacked mid-IR s pectr um of z ~ 
1.1 ULIRGs obse rved by iDasvra et all (|2009f ) with the 
local SEDs from IRieke et al.l (|2009D . Fig. g] suggests 
that the spectral features of local L(TIR)~ 10 11,3 L 
galaxies should be similar to those for z ~ 1 ULIRGs 
with L(TIR) ~ 10 12 - 3 L . We confirm this prediction in 
Fig. [7J where the average observed spectrum is indeed 
consistent with local LIRGs with L(TIR) of - 10 11 L Q 
and clearly departs from the local SED for ~ iq 12 25 L 
ULIRGs. 

These tests support our hypothesis that Xl(tir) is a 
valid tracer for the starburst environment and the ob- 
servable spectral features. More importantly, they inde- 
pendently confirm that high-redshift star-forming galax- 
ies, including SMGs, have similar star- forming environ- 
ments to local normal star-forming galaxies and their 
higher star-formation rate is primarily due to higher 
^sfr. 

3.2. A Unified View of Star- Forming Galaxies 

We believe from the nearly consistent sizes of local and 
high-redshift star-forming galaxies that the most signif- 
icant evolution between these two populations is in the 
Ssfr- In other words, while the sizes of these galaxies 
remain about the same, their star formation rate den- 
sities increase greatly from normal star-forming galaxies 
to SMGs. 

Apart from many high-resolution observations of high- 
redshift SMGs that find them to be physically extended 
systems, recent observations of the 158 fim [CII] line 
provide another probe into the environment of the star- 
forming regions. The [CII] line is an important cooling 
line for the photodissociation regions at the surfaces of 
molecular clouds. Combining this [CII] line with the CO 
(1 — 0) line yields a color-color diagram of L[ CII ]/Lir vs. 
L CO (i-o)/Lir that is sensitive to both the incident UV 
flux and the gas density. lHailev-Dunsheath et al.l (j2010h 
study the environment of the z = 1.3 galaxy MIPS 1428 
that has L(TIR) ~ 10 13 L and report that its L[c n i/LiR 
is a factor of ~ 4 higher than for local ULIRGs while the 
Lco(i-o)/Lir ratios are comparable, indicating a simi- 
lar incident UV flux in both populations but that the 
gas density of MIPS 1428 is ~ 100 x lower than those 
in local UL IRGs. This b ehavio r suggests a galaxy- wide 
starburst. Ilvison et al.l (j2010cf ) apply this analysis us- 
ing Herschel measurements, resulting in a similar find- 
ing for SMMJ2135 at z = 2.3. The star- forming envi- 
ronments in these two high-redshift galaxies appear to 
be similar to M82 and other n ormal starburst galaxies 
rather than to local UL IRGs (jHailey-Dunsheath et al.l 
120101 : Ilvison et al.ll2010cl ). SMMJ2135 is a particularly 
notew orthy case becau s e it is lensed by 32 x, which al- 
lowed iSwinbank et al.l (j2010b[ ) to observe it with the 
SMA at 0^3 x 0^2 resolution. They resolved the galaxy 
into four ~ 100— pc massive star-forming regions dis- 
tributed across the projected distance of 1.5 kpc, directly 
confirming the distributed nature of the galaxy- wide star 
formation. 

The optical morphologies of the subsample observed 
by GOODS, shown in Fig. [2j suggest that - 2/3 of the 
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subsample are quiescent, normal galaxies and the other 
~ 1/3 show signs of disturbed morphologies. The frac- 
tion of disturbed morphologies increases with redshift. 
However, only ~ 5 systems out of 36 can be identi- 
fied positively as interacting systems, while the rest of 
those with disturbed morphologies could as well be due 
to instabilities fueled by rapid, asymmetric infall of gas 
resulting in large clumps of star- forming region s simi- 
lar to those seen in SMMJ2135. iLehnert et all (|2009f ) 
also observed clumpy, galaxy- wide starbursts in 11 star- 
forming galaxies at z ~ 2 in rest-frame optical wave- 
lengths. The optical surface brightnesses for their sam- 
ple are more than an order magnitude greater than for 
local star- forming galaxies, consistent with our result. 

The molecular gas and sta r formation relation, the 
Kennicutt- Schmidt Law (e.g. Schmidt] 119591 : iKennicuttl 
1998), is shown bv iGenzel et al.l (|2010f ) to have a slope 
of 1.1 to 1.2 over a large range of stellar mass surface den- 
sity (lO 5 to 10 4 M pc~ 2 ) for both low and high-redshift 
samples. A remarkable difference between low and high- 
redshift star-forming galaxies is that the gas depletion 
time increased from 0.5 Gyr at z ~ 2 to 1.5 Gyr locally 
(jGenzel et al.1120101 ). which is consistent with the picture 
that star-forming galaxies at low and high-redshift har- 
bor similar star-forming environments but the gas con- 
sumption rate, and hence the star formation rate, is sig- 
nificantly higher at high-redshift. 

The finding that the physical conditions in high- 
redshift galaxies' star- forming regions are similar to those 
in local quiescent star-forming galaxies indicates that 
their intense star formation is unlike the transient star- 
bursting phase due to rapid infall of g as as a result of 
galaxy interaction, as seen in local ULIRGs. Rather they 
may represent an isolated evolution which could be ob- 
servable for an extended period of time. This picture is 
supported by the behavior of massive star-forming galax- 
ies a^Jngli^red^hift found in the cosmolo gical s imulations 
of lAgertz et"all (|2009f ) and IDave et~aTl (|2Q10f ). The lat- 
ter simulated populations with observational properties 
consistent with SMGs consists of isolated galaxies in the 
middle of large potential wells with larg e gas r eservoirs. 
It should also be noted that IDave et all ([2010) report a 
highly asymmetric distribution of gas density, star forma- 
tion, and velocity field in the simulated SMGs consistent 
with the disturbed morphologies observed. 

4. CONCLUSION 

We made a compilation of physical size measurements 
for 43 local galaxies with L(TIR) ranging from normal 
star-forming galaxies at 10 9 Lq to ULIRGs at > 10 12 
L , as well as 48 intermediate and high-redshift galaxies, 
including SMGs. 

Our compilation shows that (1) the physical scale of 
high-redshift ULIRGs and submillimeter galaxies is con- 
sistent within an order of magnitude with that of local 
normal star-forming galaxies (4.4 kpc in median diam- 
eter), while local LIRGs and ULIRGs are significantly 
smaller (0.8 kpc in median diameter); (2) there is a cor- 
relation of L(TIR) and Dl(tir) extending over five orders 
of magnitude in L(TIR) for normal star- forming galax- 
ies and high-z galaxies. Local LIRGs and ULIRGs have 
significantly higher El(tir) than high-redshift galaxies 
with similar L(TIR) and diverge from this correlation. 
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Fig. 6. — The logarithmic ratio of the observed 24 jim. an d 1.4 
GHz flux ratios and those predicted by the Rieke et al. (2009) SED 
templates based on local star-forming galaxies. Color coding rep- 
resents the L(TIR) for each object as in Fig. [5] The increase of 
the ratios above z ~ 1 indicates that the aromatic emissions at 
this redshift range are stronger than expected based on local SED 
templates. 
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Fig. 7. — The stacked observed spectrum of ULIRGs at 0.3 < 
z < 3.5 with aver age redshift z = 1.1 (thick blac k line) from 
IDasvra eF al. (2009) compared to the SED templates of I Rieke et alJ 
(2009) derived for local galaxies (color coded lines). This plot is 
normalized by the flux at the observed 24 fim bandpass of Spitzer 
at z = 1.1 (red dotted line). The z ~ 1 stacked ULIRG spectrum 
with L(TIR) ~ 10 12-3 L© has spectral features consistent with 
local galaxies that have L(TIR) ~ 10 11-3 L©, agreeing with the 
prediction from the relationship in Fig. 3] 

The fact that we do not find a significant devia- 
tion from this relationship in high-redshift galaxies with 
L(TIR) in the ULIRG range indicates that the local 
ULIRGs as well as LIRGs belong to a rare population 
likely driven by a unique process. High-resolution studies 
of local ULIRGs have pointed out their disturbed mor- 
phology, double nuclei, and other signs of merger activ- 
ity. The correlation we have found, however, suggests 
that the high L(TIR) of SMGs and ULIRGs at large 
redshifts can be explained, to first order, by the higher 
^sfr within isolated, quiescent galaxies. 



MORPHOLOGY AND SIZE OF LUMINOUS INFRARED GALAXIES 



9 



We thank Anita Richards and Tom Muxlow for ra- iMenendez-Delmestre et all ((2009). We also thank Ben- 
dio data in the HDF and insig htful discussions ; Scot t jamin Weiner for his assistance with data compilation 
Chapman for radio images from lChapman et al.l (12004); and valuable discussions. This work is supported by con- 
Kalliopi Dasyra for spectra from Das vra et al.l (j20 09); tract 1255094 from Caltech/JPL to the University of Ari- 
Karin Menendez-Delmestre for spectra published in zona. WR gratefully acknowledges the support from the 

Thai Government Scholarship. 

APPENDIX 

A. ESTIMATION OF PHYSICAL SIZES FOR LOCAL IR GALAXIES 

Our local sample of star-forming galaxies was selected to represent a broad range of IR luminosities from normal 
star-forming galaxy (L(TIR) < 10 11 L©), luminous IR galaxy (LIRG, 10 11 < L(TIR) < 10 12 L ), to ultraluminous IR 
galaxy (ULIRG, L(TIR) > 10 12 £©). Here we discuss how we select the local starburst galaxies and the methods we 
employ to estimate the physical sizes to compare with high-redshift star-forming galaxies in the same metric. For all 
but t hree local galaxies, L(TIR) was obtained from the IRAS Revised Bright Galaxy Sample (RBGS; Sanders et al.l 
I2003D . The RBGS provides L(TIR) that has a definition similar to that of L(TIR) from the iRieke eFall ((2009) 
formalism used for our intermediate and high redshift compilation. The three exceptions are discussed separately. 

A.l. Normal Star- Forming Galaxies 

Galaxies from the Spitzer Infrared Nearby Galaxies Survey (e.g., iKennicutt et al.l [20031 : iCalzetti et al.l [20071 ) are 
used to represent normal star-forming galaxies. SINGS galaxies are selected from within the local volume to allow IR 
imaging at reasonably good physical resolution. We use the 24 /am images of these galaxies to map star formation. 
We only select high-luminosity galaxies. As noted by ICalzetti et al.l ([20071 ), NGC 4125 and NGC 5195 contain Seyfert 
2 nuclei and thus are excluded from our compilation. L(TIR) for the other galaxies was taken from the RBGS except 
for NGC 1512, NGC 2841, and NGC 4625 that are not in the ca talog. For these th ree galaxies we use 24 /im fluxes 
from lDale et al.l (|2007D to estimate L(TIR) using the appropriate IRieke et all ([2009) SED. 

We convolved the SINGS 24 /im MIPS images with a series of Gaussians that have angular FWHM corresponding 
to physical sizes ranging over 1 — 8 kpc at the galaxy's distance and then fit 2D Gaussians to estimate a deconvolved 
size from the convolved images. We found that the recovered FWHM is a slowly varying function of the convolved 
Gaussian's FWHM, which assures that sizes estimated from this method are robust. For the actual size measurement 
for these galaxies we convolve all SINGS 24 /im imaging with a Gaussian that has FWHM corresponding to 4 kpc 
at each galaxy's distance and measure deconvolved sizes from the convolved image. This procedure is illustrated by 
NGC 3627 in Fig. [TJ These sizes are reported in Table [1] for normal star-forming galaxies. 

We have 19 galaxies from the SINGS sample with L(TIR) in the range of 10 8 6 - 10 10 6 L©, a median L(TIR) of 
10 9 ' 8 L , and a median physical diameter of 3.5 kpc. 

A.2. LIRGs 

Our local LIRGs were drawn from lCondon et al.l (|1991f ). lAlonso-Herrero et al.l (|2006l ). and llono etall (|2009l ). Their 
size measurements are based on the 8.4 GHz VLA radio observations, Paa imaging using Hubble Space Telescope 
( HST), and Submillim eter Array (SMA) submillimeter observations, respectively. 

iCondon et al.l (|1991f ) observed 40 LIRGs selected from the IRAS Bright Galaxy Sample using the VLA with a 
resolution of 0'/25. From th e se 40 LIRGs, 15 are dominated by a compact radio component likely to be an AGN, 
referred to by ICondon et al.l ([19911 ) as "monsters" and another 5 are known to harbor AGN according to SIMBAD, 
which we exclude; 6 are multi-component with one or more components not reported, which we exclude; 2 are multi- 
component with well c onstrained sizes f or both components (NGC 3690, which is Arp 299 for which we adopt CO 
size measurement fro mllono et al.l (120091) . and IRAS F15163+4255) ; one has an obse rvational issue (NGC 1614) but it 
was also observed by Ion o et all ([2009) and therefore we adopt the llono et al.l ([2009) size measurement for NGC 1614; 
there are size measurements for only four galaxies of the remaining 11. These four are UGC 2369, IRAS F03359+1523, 
UGC 4881, and IRAS F17132+5313. Along with the aforementioned I RAS F15163+4255, a multi-component object 
with well -measured size, we have f ive LIRGs from ICondon et all (fl99ll) . 

The lAlonso-Herrero et al.l (|2006[ ) LIRGs were selected such that their Paa emission line would fall in the narrow band 
of the F190N filter of the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) on HST. The resolution 
of NICMOS (O'.'076/pixel) resolv ed detailed structures of these LIRGs and thus we need to measure their sizes using 
the same procedure as in § IA.1I However, we convolved the images using Gaussians with FWHMs corresponding to 
only 1 kpc at each LIRG's distance. This was necessary because of the intrinsically smaller sizes of the LIRGs, which 
we also tested by measuring the diameter encircling 90% of flux. The results from our measurements are given in the 
L IRGs sectio n of Table [J 

llono et al.l (2009) report sizes of LIRGs within 200 Mpc using the SMA, based on CO (3 — 2) observations by 



Wil son et al.l (|2008[ ). The sample contains LIRGs harboring AGN (including NGC 2623, NGC 6240, UGC 5101, Mrk 
231, Mrk 273, IRAS F10565+2448 and IRAS F17207-0014) and multicomponent objects where the sizes of some 
components cannot be determined (including Arp 55, Arp 299, NGC 5331, and NGC 5257/8) and thus we are left 
with four starburst-dominated objects with good size estimates. They are Arp 193, VV 114, NGC 1614, and IRAS 
10565. 
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In total we have 21 LIRGs with median L(TIR) of 10 11 * 5 Lq and a median physical diameter of 0.8 kpc. This 
sample is dominated by the 12 galaxies with Paa images, and these galaxies by themselves should provide an unbiased 
estimate of sizes of the star forming regions in LIRGs. The estimated sizes of the galaxies measured in the radio and 
in CO are consistent with those from Pa<x 

A.3. ULIRGs 

Despite the many local ULIRGs, selecting a sample to study the physical sizes of the starburst regions is complicated 
by two issues. First, a majority of these ULIRGs harbor AGN and thus do not represent a starburst environment. 
Second, the sizes of the IR-emitting regions are significantly smaller than the optical and near-IR extents of these 
galaxies and the existing mid-IR data in the literature (e.g. SINGS) do not resolve the ULIRGs' IR emission. Our 
compilation therefore relies on radio and submillimeter interferometric observations of four ULIRGs that are known 
to have their IR luminosity dominated by star-fo rming emission : VII Zw 31, IRAS F23365+3604, Arp 220, and IRAS 
F17207-0014. Their IR luminosities according to iSanders et al.l (j2QQ3h are in the range of 10 12 00 - 10 12 45 L . 

Arp 220 has two nuclei ~370 pc apart; the western nucleus is slightly more luminous than the eas tern nucleus 
(lR,ovilos et al.ll2003[ ). This galaxy was studied with Very Long Baseline Interferometry (VLBI) at 18 cm bv lSmith et al.l 
(1998) who found the western nucleus to contain most of the individual 18 cm sources (13 sources in the western nucleus 
vs. three in the eastern nucleus). They report that the sources in the western nucleus are confined within 75 pc xl50 
pc. Circularizin g these value s give s a diameter of 0.1 kpc for the western nucleus alone. The 1.6 and 5 GHz MERLIN 
radio maps by iRovilos et al.l (j2QQ3h indicate both nuclei have similar physical sizes. We estimate the physical size of 
Arp 220 from the 5 GHz map where the two nuclei are clearly separated. Combining the emission regions from both 
nuclei, we adopt an effective diameter of 0.2 kpc for the starburst r egion of Arp 220. 

The sizes of VII Zw 31 and IRAS F23365 +3604 are r eport ed bv iDownes fc Solomon! (|1998[ ) based on CO (1 - 0) 
observations; the IRAS F17207-0014 size bv llono et al.l ((2009) is based on CO (3 — 2) observations. Again, we treat 
these CO-derived sizes as upper limits for the size of the IR-emitting extent of the galaxy. 
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TABLE 1 

Local Compilation of Star-Forming Galaxies 



Source 


IRAS ID 


Distance a 


L(TIR) a 


Diameter References 13 






(Mpc) 


OgL©) 


(kpc) 



Normal SF Galaxies 



NGC 2976 


F09431+6809 


3.8 


8.61 


1.5 


1 


NGC 4826 


F12542+2157 


6.0 


9.14 


1.2 


1 


NGC 2403 




Q Q 

o.o 


y.zo 


o.U 


i 


NGC 925 


F02242+3321 


9.8 


9.46 


6.1 


1 


NGC 


1512 




11.3 


9.49 


3.2 


1 


NGC 


5866 


F15051+5557 


13.0 


9.51 


1.8 


1 


NGC 


2841 




10.5 


9.52 


4.9 


1 


NGC 


4559 


F12334+2814 


11.9 


9.62 


5.0 


1 


NGC 


4736 


12485+4123 


5.7 


9.79 


1.7 


1 


NGC 


3198 


F10168+4547 


14.7 


9.81 


3.8 


1 


NGC 


3184 


10152+4140 


11.9 


9.86 


6.5 


1 


NGC 


3351 


F10413+1157 


10.8 


9.89 


2.2 


1 


NGC 


3938 


F11502+4423 


13.1 


9.99 


4.7 


1 


NGC 


4569 


F12343+1326 


17.8 


10.08 


3.1 


1 


NGC 


5055 


F13135+4217 


8.4 


10.15 


4.3 


1 


NGC 


5033 


F13111+3651 


14.7 


10.19 


1.1 


1 


NGC 


3627 


F11176+1315 


9.3 


10.44 


3.5 


1 


NGC 


5194 


F13277+4727 


8.8 


10.48 


5.7 


1 


NGC 


7331 


F22347+3409 


16.2 


10.64 


4.7 


1 


LIRGs 














NGC 23 


F00073+2538 


63.9 


11.11 


1.2 


2 


NGC 6701 


F18425+6036 


60.6 


11.11 


0.7 


2 


UGC 1845 


F02208+4744 


66.4 


11.13 


0.81 


2 


NGC 5936 


F15276+1309 


65.1 


11.13 


0.56 


2 


MCG +02-20-003 




TO A 
I Z.4 


ii i/i 

11. 14 


U. ( o 


Z 


NGC 2369 


F07160-6215 


47.1 


11.16 


0.82 


2 


ESO 320-G030 


F11506-3851 


40.4 


11.16 


0.91 


2 


IC 5179 


F22132-3705 


50.0 


11.22 


1.6 


2 


NGC 2388 


F07256+3355 


61.9 


11.29 


0.84 


2 


NGC 7771 


F23488+1949 


61.2 


11.40 


1.04 


2 


MCG +12-02-001 


F00506+7248 


68.9 


11.50 


0.8 


2 






F03359+1523 


146.9 


11.53 


0.07 


3 


NGC 1614 


F04315-0840 


67.1 


11.66 


0.9 


4 


UGC 2369 


F02512+1446 


130.7 


11.66 


0.08 


3 


Arp 236 


F01053-1746 


84.2 


11.71 


2.6 


4 


Arp 193 


F13182+3424 


107.1 


11.73 


0.8 


4 


UGC 4881 


F09126+4432 


172.7 


11.75 


0.14 


3 


Arp 299 


F11257+5850 


51.2 


11.94 


0.42 


4 






F17132+5313 
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a Distance and L(TIR) from Sanders et al. (2003) and adjust to match our cosmology. 
b Diam eters references 1. m e asured from Spitzer MIPS 24 fim imaging taken by 
SINGS (|Kennicutt et al.l 2003; Cal zetti et all 120071, e.g.,), 2. measured from Hubble 
NICMOS Pa-a taken by [Alonso-Herrero et al. (2006), 3. 8.4 GHz rad io sizes given by 
ICondon et all (1991), 4. C O (3 - 2) sizes given bvllono et al.l (|2009), 5. CO (2 - 1) 
or CO (1 — 0) sizes given bv lDownes fc Solomon! (1 19981 ), 6. 5 GHz radio size based on 
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Daddi et al. (2010), 10. GOODS Spitzer Legacy Data, Dickinson et al in prep., 11. Ilvison et al.l (|20071 ). 12. 
Muxlow et al. (2005), 13. Weiner, private comm., 14. IMorrison et al.l (2010) 
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